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Abstract 

We start with the total energy E for a system of three scalar relativistic 
particles that, because of Einstein’s relation, will have square roots of functions 
of the momenta. By taking powers of this relation, we finally get a fourth degree 
polynomial in where the square roots have disappeared, and which we can 
convert into a type of Schroedinger equation. To be in the center of mass frame 
we pass to Jacobi momenta and then replace them by creation and annihilation 
operators. We thus get an equation in terms of the generators of a U(2) group, 
which, in principle, we can solve in an elementary way. Finally we rewrite our 
equation in a Poincare invariant form. 


1 Introduction 

In the II Harmonic Oscillator Conference I presented a paper dealing with systems 
of relativistic particles interacting through Dirac oscillators. The results were later 
applied to the mass spectra of baryons and mesons [1,2,3]. 

As all the results presented had already been published I prefer to deal in this paper 
with a new approach, restricted here to scalar particles, that seems to me a systematic 
way to attack many body problems with oscillator interactions. 
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I will start by considering a non relativistic problem of n free particles and indicate 
the steps by which later it can be reduced to a system with oscillator interactions, 
which will serve as a model for the relativistic problem we wish to analyze. 


2 The non-relativistic problem 


Let us first start with a system of n free non-relativistic particles of the same mass m, 
and take units 

= m = c = 1, (2-1), 

where the velocity of light will appear only in the next section, but we want to start 
from the beginning with units in which everything is dimensionless. 

The classical total energy is then 

« = (1/2)Ep.-P. (2-2) 

where p, are the three dimensional classical momentum vectors of particle s. 

From the beginning we would like to work in the center of mass frame, because 
our interest will be the internal energy of the system and not the contribution from its 
center of mass motion. The best way to achieve this is to pass to Jacobi coordinates [4] 
defined by the orthogonal transformation 


P', = [«(« + !)] 2 ^pt-sp,+i],s = l,2,...n- 1, 
_i * 

Pn = « 2 ^p„ 


(2.3a) 


(2.36) 


Clearly p(, is proportional to the total momentum and in the center of mass system 
it will vanish, so Eq. (2.2) reduces to 




(2.4) 


»=i 


The Schroedinger equations corresp>onding to (2.4) is obtained when we replace p' 
by the operator 


n' = iA 

*dx;’ 


(2.5) 
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with being the corresponding Jacobi coordinate vector. As the are hermitian 
operators we can also write Schroedinger equation as 

lY. (p's ' p's)^ = 

5=1 ^ ^ 

We can easily transform this into an hermitian oscillator operator equation if we 
make the replacement 

p; p; - tu;x;, (2.7a) 

p's p's + ( 2 - 7 ^) 

where the second equation follows from the first as both Pj,x^ are hermitian operators. 

Thus we now get a Schroedinger equation of the form 

I E \p's + - (3/2)o;(n - 1)1 0 = (2.8) 

S=1 

whose eigenvalue for the energy E will be 

E = uN, (2.9) 

with N being the total number of quanta i.e. 

( 2 . 10 ) 

5=1 

The previous analysis is standard except for the fact that we start from a system 
of n free particles. Furthermore our notation in terms of three vectors and Jacobi 
coordinates, avoids the worry about the Galilean invariance of the whole procedure. 
We will now consider a similar set of steps for a relativistic problem. 

3 The system of three relativistic particles 

Rather than discuss the system of n relativistic particles, we shall restrict ourselves to 
n = 3, as we will see that the case is general enough, with only the algebraic steps 
becoming more complicated as n increases. 

In our units the total energy for a system of three free relativistic particles can be 

E = ±Ux±U2± Ra, 
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written as 


(3.1) 



where = 1,2,3 is defined as 


n. = (p; + i)l. (3.2) 

It is very important to note that, in our units, Einstein relation is + 1, and 

when reduced to the E itself gives both the square root in (3.2) and the ± signs in 
(3.1). 

Obviously we can not get a Schroedinger equation from the relation (3.1), but we 
can take ills to the right hand side and square both sides. Then we can square again 
and again appropriately, and we easily arrive at the fact that (3.1) becomes an eight 
degree equation in E (actually of fourth degree in E^) of the form 

^E\n]) = 

E^ - AAE^ + (4/1" + 2B)EA - (4(7" + AAB)E'^ + 5" = 0 (3.3) 

where A, B, C are functions of H", s = 1, 2, 3 given by 

A = Ul + Ul + II3, (3.4a) 

5 = - 2n^n^ - 2n^n^ - 2njn^, (3.46) 

(7" = lenjn^ri^. (3.4c) 

Now we can write an equation that does not have E as an eigenvalue, but in which 
it appears as a parameter, if we replace by —id/dxs as in (2.5), so that 11^ become 
the operators 

= + (3.5) 

and we get 

$(E", n,")V’ = 0 (3.6) 

Thus far we have obtained nothing useful because Pi, P2, Ps considered as operators 
of the form (2.5), commute with the operator $ and so are integrals of motion, so that 
Ip can be written as 

Ip = exp[i(pi • Xi + P2 • X2 + P3 • X3)], (3.7) 

where now pi , P2, Ps are ordinary numbers and we are returned to equation (3.3) whose 
eight roots for the energy E are obviously given by (3.1) with all the possible combi- 
nation of the signs ±. 
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Before proceeding further, along the lines of the previous section, we again remark 
that we would like to work in the center of mass frame, as our interest is restricted to 
the internal energy of the system. Thus we go, as in section 2, to the Jacobi momenta 
p',s = 1,2,3 which from (2.3) are given now by the matrix relation 

/P'l \ / ^ ® \ /Pi \ 

JL J_ _ /^ 

V 3 

J_ J_ J_ / 

>/3 \/3 ' 

As the matrix is orthogonal, transposing it we get p, in terms of p' , and as we want 
to be in the center of mass frame pg = 0, so we get finally 


P2 


P2 

\P3/ 


(3.8) 




1 


1 


^P'l ^P2.P3 


^P2> 


(3.9) 


As Eq. (3.6) contains only powers of 11^,5 = 1,2,3, we can write the latter using 
the hermitian property of p', now considered as operators of the type (2.6), as 

n? = ^P? • pi + ^P2 • ?2 + ^(Pi^ ■ P2 + P2^ ■ Pi) + (3.10 g) 

A 2 = \pi • pi + ^p' 2 ^ • p '2 - ■ P 2 + P 2 ■ pi) + (3.106) 

n^ = ^p'2'p'2 + i (3- 10c) 

The interesting point is to introduce the oscillator interaction, exactly as in the 
replacement we made in (2.7) in the non- relativistic problem. For notational purposes 
we introduce the creation and annihilation operators 

11 1 

»?. = - iu; 2p^),^ = l,2, (3.11a) 

C = ;^(^2xi + iu;"2p'J,5 = 1,2, (3.116) 

so that the relations (2.7) can be written as 

pi -iu;2v/2C, (3.12a) 

p9 -> ?w2v^i;,, (3.126) 

Under this replacement the operators become then 

fli = c<;[Cn + 3^22 + ^ C 21 )] + 1 
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(3.13a) 



(3.136) 


nl = + Cj,)! + 1 

^3 = -^(^C 22 + Ij (3.13c) 

where the operator Cst, s,t = 1,2 are defined by 

Cat — Va ' ( 3 - 14 ) 

From the fact that 

~ ~ 1»2,3 J S = 1,2, (3.15) 

we have the commutation relations 

[Cat-, Ca't>] = Cat'^a't ~ Ca't^at'i (3.16) 

and thus they are generators [5] of a U(2) group. Therefore the operators nj,fl 2 ,n|, 
appearing in the equation (3.6), are linear functions of the generators of this group. 


To obtain from Eq.(3.6) the eigenvalues of the energy for this relativistic oscillator 
problem we can proceed as follows: First we note that the first order Casimir operator 
of U(2) group is 


N = CnFC22, 

(3.17) 

ajid that it has an SU(2) subgroup whose generators are 


F+ = Ct2, 

(3.18a) 

> 

III 

N-» 

1 

(3.186) 

F-=C2U 

(3.18c) 

with a corresponding Casimir operator of the form 


= F_F+ + Fo{Fo + 1). 

(3.19) 


A A 

The by definition commute with all Cat and among themselves, so from (3.13), 
they will be integrals of motion of the operator $(£^,11^). Thus the eigenstates of the 
Eq. (3.6) can be characterized by the eigenvalues of TV, which we denote respectively 

by 




(3.20) 
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with / taking the values (iV/2, (A^/2) - 1, down to 1/2 or 0 depending on whether N 
is odd or even. 

Another operator that commutes with N, F'^ is obviously and we shall designate 
its eigenvalue by 

IX = /, /-I,. (3.21) 
so the eigenstates associated with N ^ F^, Fo could be represented by the ket 

\Nfu >, (3.22) 


and the solution V’ of Eq. (3.6) is necessarily a linear combination of these kets i.e. 

^ ^ a,\Nfu >, (3.23) 

as N , F^ are integrals of motion. 

To obtain the eigenvalues of the internal energy E as function of N, f we need first 
to consider the matrix elements of the operator ^ of (3.6) in the basis (3.22) i.e. 

< iV/l^'|$(F^^^)|A/t/>, (3.24) 


which from (3.13), (3.17), (3.18) is a straightforward, but laborious, calculation of the 
type familiar in angular momenturh theory, as the group there is also SU(2). 

To get the internal energy 

E(NJ,a), (3.25) 

with a indicating the rest of the indices, we need to evaluate the determinant of the 
(2/ + 1) X (2/ + 1) matrix whose elements are (3.24) and equate it to zero. This gives 
us a numerical equation of degree 4(2/ + 1) in the variable E^ and its solution provides 
us with values indicated symbolically in (3.25). 

As our purpose is to provide the method of solution for the internal energy of rel- 
ativistic three body oscillator problem, we will only carfy the calculation of (3.24) for 
the single case when 

A = / = 1/ = 0 (3.26) 

which imphes that 

< 000|n,^|000 >= 1, s = 1,2,3 (3.27) 

and so A, B, in (3.4) become respectively 


A = 3, 


B = -3, 
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= 16, 


(3.28) 



and the equation for the energy is given by 

- 12E^ + - 2^E^ + 9 = 0, (3.29) 

whose four roots for E^ are E^ = \^ repeated three times, and E^ = 9, with E = ±3 
and ±1 as we expect from (3.1). 

So far we have discussed, and given a method for solving, the equation related with 
a three body relativistic problem with an oscillator interaction. In the next section we 
proceed to show that we can formulate it in a Poincare invariant form. 


4 Poincare invariance of the three body relativis- 
tic equation with oscillator interactions 

To express Eq. (3.6) in a Poincare invariant form we start with definition of the total 
four momentum for the three particle problem i.e. 

Pt, = PmI + Pm2 + (4. 1 ) 

where // = 0, 1,2,3 with Posi^ = 1,2,3, being the time like component while = 
1,2,3, are the space like components of the vector p* of the previous section. 

We shall require also a unit time like four vector which we shall define as 

= Pui-g'"'^ P<tPt)~^ ^ (4.2) 

where repeated indices <t, r are summed over 0,1,2, 3 and our metric tensor is taken as 
= 0 if <T / T, gu = 922 =933 = -9oo = 1 (4.3) 

Clearly in the center of mass frame where P, = 0, r = 1,2, 3, takes the value 

(u,) = (l,0,0,0) (4.4) 

The operators — 1,2, defined in (3.11) are space like three component 

vectors which could be denoted by r}ia,Cis,i = 1,2,3. A time like component could be 
added through the definition (3.11) just by putting p'^stPos instead of and thus 

we would get Posiios which, together with rfis,iis, form the four vectors 


P 6 , 1 , 2 , 3 , 5 — 1 , 2 . 
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(4.5) 



We do not want. to use these operators directly in the definition of the of (3.14), 
but rather utilize their transversal parts defined by 

its = + {g^^i<rsUr)u,. (4.66) 

These transverse operators have the property that in the center of mass frame where 
= (1000) we have, because of the matrix (4.3), that 

nts = 0 . its = 0 . gts = ms. its = ( 4 - 7 ) 

Thus now the generator C,t,s,i = 1,2, appearing in the definitions (3.13) of 
nf , ftj, n| can be expressed in a Lorentz invariant way by 

Cst = g^'^vtsitv (4.8) 

as in the center of mass frame it takes the form (3.14) i.e. Cst — Vs ’C- 

As for the energy appearing in Eq. (3.6) it can be substituted by the operator 

P,Pr). (4.9) 

because in the center of mass frame F, = 0, i = 1, 2, 3, and from the metric tensor (4.3), 
we see that the parenthesis in (4.9) reduces to Fj, which is the time like component 
of the four momentum vector squared and thus corresponds to the square of the total 
energy of the system. 

With the definitions (4.8) of Cst and (4.9) of E'^ substituted in Eq.(3.6) we get a 
Poincare invariant equktion for our problem, as Cst. given in terms of Jacobi cejordinates 
and momenta, is also invariant under translation in space time, and thus commutes with 

P.- 

We have then arrived at a procedure for deriving a Poincare invariant equation for a 
three particle system with oscillator interactions which, in the center of mass reference 
frame, can be solved by a simple group theoretical procedure, which leads eventually 
to algebraic equations of degree 4(2/ "hi) for that can be solved numerically to 
give the spectrum of the problem. 
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